The Spatially Resolved Ha-Emitting Wind Structure of P Cygni 

Aurelian Balan^'^, C. Tycner^, R. T. Zavala^, J. A. Benson^, D. J. Hutter^, M. Templeton^ 

ABSTRACT 

High spatial resolution observations of the Ha-emitting wind structure associ- 
ated with the Luminous Blue Variable star P Cygni were obtained with the Navy 
Prototype Optical Interferometer (NPOI). These observations represent the most 
comprehensive interferometric data set on P Cyg to date. We demonstrate how 
the apparent size of the Ha-emitting region of the wind structure of P Cyg com- 
pares between the 2005, 2007 and 2008 observing seasons and how this relates to 
the Ha line spectroscopy. Using the data sets from 2005, 2007 and 2008 observ- 
ing seasons, we fit a circularly symmetric Gaussian model to the interferometric 
signature from the Ha-emitting wind structure of P Cyg. Based on our results 
we conclude that the radial extent of the Ha-emitting wind structure around 
P Cyg is stable at the 10% level. We also show how the radial distribution of 
the Ha flux from the wind structure deviates from a Gaussian shape, whereas 
a two-component Gaussian model is sufficient to fully describe the Ha-emitting 
region around P Cyg. 

Subject headings: stars: winds, outflows — stars: individual (P Cyg) — tech- 
niques: interferometric 

1. Introduction 

Luminous Blue Variable (LBV) star P Cyg (HD 193237, B2pe) is an unusual star with 
a unique stellar wind structure. LBVs are evolved, very luminous, massive supergiant stars 
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t hat show some type of instabihty. As the hkely progenitors of Wolf-Rayet stars (ICrowther 
20071 ) LBV's may provide insight for the ultimate end-stage of these stars, the core-collapse 
supernovae. P Cyg was discovered ii i 1600 after a vi olent mass- loss event caused it to quickly 
brighten to a third magnitude star (jde Grootlll969l ). The violent mass-loss events that spur 
these eruptions are not fully understood, but a number of cor npeting models including singl e 
and binary star scenarios attempt to explain these eruptions ( Humphrevs fc Davidson 1994[). 



With a mass-loss between 2x10 ^ and 4x10 ^ Mq yr ^ (Ivan Blerkom 



1978 



Abbott et al. 



1980l : iLeitherer fc Zickgrail 119871 ) . P Cyg displays characteristics that are suggestive of en- 
ergetic mass outflows. It is for this reason that P Cyg has a unique stellar wind structure 
with possible high density regions in motion within the wind structure. Despite numerous 
studies published in the literature that discuss P Cyg or other LBVs, there is still need for 
observations that can spatially resolve the circumstellar region around P Cyg. 

The spectrum of P Cyg reveals emission line profiles with their archetypical shapes, 
which are formed by the circumstellar matter that has been shed from the central star. These 
profiles are commonly identified by an emission component to the red side of an absorption 
line. The origin of this Doppler shift can be explained with a sphericall y symmetric outflow- 
ing w ind in which the velocity increases with radial distance (see, e.g., iLamers fc CassineUi 
19991 . § 2.2). The absorption component corresponds to the continuum light absorbed by the 
wind structure directly between the observer and the central star. The emission component 
on the other hand originates in the spherically symmetric halo around the central star. Be- 
cause the same characteristics apply to the Ha emission line, which is one of the strongest 
emission lines in the spectrum of P Cyg, this makes the Ha emission line an excellent probe 
of the wind region around P Cyg. 



P Cyg currently has an estimated mass of 30 ± 10 Mq (ILamers et al. 



1983 



1985). At the 



start of its lifetime P Cyg possibly had an initial mass of 50 Mq (jde Groot fc Lamerslll992l ). 
but has shed much of its mass because of its violent history of mass loss events. The effective 
temperature and the stellar diameter of P Cyg a re estimated to be Teg = 19,300 ± 2000 K and 
76 ± 14 i?0, respectively (ILamers et al.lll983l ). However, the Ha-emitting region extends 
radially much farther than the central star allowing us to resolve this region using long- 
baseline interferometric techniques, even though the distance is estimated to be in the range 
of 1.7 ± 0.1 kpc to 1.8 ± 0.1 kpc jLamers et al.lll983l : iNaiarro et aDll997h . 



The circumstellar region of P Cyg has been spatially resolved in the past using ra- 
dio and optical interferometry, as well as direct imaging with adaptive optics (AO). Radio 
interferometric observations detec t the nebula aroun d P C yg at the angular scales from 
~50 mas to almost an arcminute (jSkinner et al.l 119971 . Il998l ). whereas optical long-baseline 
interferometry is sensitive to structures in the nebula of P Cyg that are more than an or- 
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der of magnitude smaller. Therefore, optical interferometry provides a unique window of 
opportunity to resolve the inner portion of the wind structure. 

Most relevant to our study are the results that probed the Ha-emitting region of P Cyg. 
For example, bas ed on near-diffra c tion l imited observations obtained using AO system on a 
1.52 m telescope, IChesneau et al.l ( 120001 ) not only resolved the outer Ha-emitting region of 
the extended envelope, but detected signatures of clumping. Although, the angular scales 
sampled with a 1.52 m telescope were quite large, in the region of 600 mas, the a ngular 
resolution was limited to 50 mas. On much smaller angular scales, IVakili et al.l ( 119971 ) using 
a single 17.7 m baseline of the GI2T interferometer have resolved the inner Ha-emitting 
envelope of P Cyg yielding an angular size of 5.52 ± 0.47 mas and reported a marginal 
detection at the Hel 667 8 line. The spatial scales probed by our study are similar to those 
sampled by IVakili et al.l ( 119971 ). however, the interferometric observations at the Ha line 
obtained in our study were acquired using 11 unique interferometric baselines ranging in 
length from 15.8 to 79.3 m, thus resulting in much higher angular resolution measurements. 



2. Observations 

2.1. Interferometry 

High angular resolution long-baseline interferometric observations were obtained on 
P Cyg over thirty-seven nights in 2005, 2007 and 2008 using the Navy Prototype Optical 
Interferometer (NPOI) located near Flagstaff, Arizona. The NPOI is a long-baseline optical 
interferometer which consists of six 50 cm siderostat mirrors that send 12.5 cm diameter 
light beams down the optical feed system consisting of 15 cm diamet er pipes to an optics la b 



where the interference (i.e., a fringe) between the beams is recorded (lArmstrong et al.lll998l ). 
Observations covering a wide spectral band can be taken over multiple baselines (currently 
ranging from 19 to 80 m in length), and are geometrically compensated and modulated using 
vacuum delay-lines to guarantee the entire band is being observed simultaneously. 

The observing process involves measuring the fringe contrast over 15 spectral channels 
every 2 ms. Each 2 ms frame is then used to obtain a squared visibility {V^) measure (i.e., 
a value measures the fringe contrast). Aft erward, these values are incoherently averaged 



over 1 s intervals to create 1 s data points (IHummel et al.ll2003l ). The final observations 



consist of 30 s of integration (also known as a "scan") during which the squared visibilities 
are measured simultaneously at three baselines per output from the beam-combiner (with 
two outputs being used simultaneously) and 15 spectral channels covering the wavelengths 
between 560 and 870 nm. The typical observational sequence consists of a pair of coherent 
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and incoherent (off the fringe) scans on a target followed by a pair of scans on a calibrator. 
The incoherent scans are used to estimate the additive bias terms affecting the squared 
visibility measures. The process is repeated during a night for as long as the fringe packet 
can be detected on the target star. Typically, the observations stop when the optical path 
difference between the different siderostats cannot be compensated anymore. 

The Ha-emitting sources form a specific class of targets for NPOI since the entire emis- 
sion line ends up only in one out of the 15 spectral channels. This characteristic allows one 
to calibrate the V"^ values from the channel containing the Ha emission using the values 



from the neighboring continuum channels ( iTycner et al.l 120031 ). Although during this cali- 



bration process one needs to assume that the diameter of the central star at the continuum 
wavelengths is known, the final results based on the Ha emission are only weakly dependent 
on the adopted value for the angular diameter of the central star (we explore this in more 
detail in § 13. 2p . For the purpose of the calibration we adopted a uniform disk diameter of 
0.2 mas to represent the central star at the continuum wavelengths. We obtain this diameter 
based on a distance estimate of 1.8 ± 0.1 kpc and a central star radius of 76 ± 14 Rq 



( iLamers et al.l 119831 ). 



Applying the procedures described above to the entire observational data set we acquired 
for P Cyg, yielded a total of 1,534 Ha squared visibilities on 11 unique basehnes (out of a 
maximum of 15 unique baselines for a 6-station configuration). Table [1] lists the number 
of unique Ha V"^ measures we obtained for each of the nights in our observing campaign. 
Although the observational setup for each of the three observing seasons was slightly different, 
and therefore a range of different baselines were sampled each year, in total we have obtained 
observations at baselines ranging from the shortest (15.8 m) to the longest (79.3 m) possible 
at the time of our observations. More specifically, the 2005 and 2008 seasons contained 
baselines between 18.9 and 53.2 m, and the 2007 season contained baselines from 15.8 to 
79.3 m in length. Figure [T] shows the sampling of the (M,f)-plane in the Ha channel for 
P Cyg based on all the available nights in 2005, 2007 and 2008. The "arcs" seen in the plot 
are a result of the changing projection of the baseline vector on the sky during the night due 
to Earth's rotation. Three seasons of squared visibility data sampled across the (m, f )-plane 
are shown in Figure [2] plotted as a function of radial spatial frequency (a/m^ + v^). The 
calibrated Ha V"^ measures are listed in Table |21 

In addition to the V"^ measures we have also obtained closure phase quantities (sum of 
complex phases on three baselines forming a triangle) on four unique baseline triangles. By 
fitting a quadratic function to the closure phases from the continuum channels (at each scan 
and closure phase configuration separately), and then removing this quadratic trend from all 
channels, including the Ha channel, we were able to obtain the differential closure phases for 
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the Ha-emitting region. Because we expect the central star to be mostly unresolved at the 
continuum channels and if we assume that the continuum source is point-symmetric, then 
we expect the closure phases to cluster around a zero value. This allows us to inspect the Ha 
closure phases with respect to the mean zero closure phase. In order to reduce the random 
noise present in data associated with a single scan we calculated the weighted mean closure 
phases for each season and each unique configuration. In the 2007 season one of the baselines 
(AE-AN) was sampled at two independent outputs from the beam combiner and this allowed 
us to calculate two closure phase quantities containing this baseline either from one output 
beam or the other. Because the observations for the other baselines were identical, we 
decided not to average these quantities, and instead treat them separately. Figure [3] shows 
the weighted mean closure phases for all the configurations with the quadratic trend based 
on the continuum channels already removed. 



2.2. Spectroscopy 

The Ha emission line, as detected by the NPOI, falls on a spectral channel that has 
a spectral width of 15 nm. Therefore, the interferometric data obtained on P Cyg does 
not provide sufficient spectroscopic information to allow us to determine the strength or the 
shape of the Ha emission line. For this reason, we obtained complementary spectroscopy 
of P Cyg using the Lowell Observatory's John S. Hall 1.1m telescope located at the same 
observing site as the NPOI. The telescope was outf itted with the Solar-Stel lar Spectrograph 



(SSS), which is a fiber-fed echelle spectrograph ( iHall fc LockwoodI Il995l ). Spectroscopic 
data on P Cyg were obtained during the period from April 2005 to December 2007, with the 
individual nights listed in Table [31 The SSS instrument produces spectra in the Ha region 
with a resolving power of ~ 10, 000. A typical Ha line profile normalized with respect to the 
continuum level of P Cyg is shown in Figure |H The equivalent width (EW) measurements 
based on the individual spectra are also plotted in the upper panel of Figure |5l 



3. Analysis 

3.1. Interferometric Results 

There are two signatures in the squared visibilities from the spectral channel containing 
the Ha emission line shown in Figure [2J One is from the central star and the other is from 
the Ha-emitting circumstellar region. To model this data, we use the same two component 
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model as was used by iTycner et al.l (120061 ) of the form: 



model 



[CpVp + (1 - Cp)Knv]^ 



(1) 



where Vp is the visibihty function representing the photosphere of the central star, Venv 
is the visibility function representing the circumstellar envelope, and Cp is the fractional 
contribution from the stellar photosphere to the total flux in the Ha channel. 

Based on the results from differential closure phases (recall § I2.ip . which showed very 
weak (if any) variations across the Ha channel, we expect the Ha-emitting structure to be 
sufficiently well described (at the level of angular resolution provided by our range of base- 
lines) by a model that is both point symmetric and concentric with the central star. This 
means that Vp and Venv in equation [1] can be represented by real functions (i.e., the Fourier 
transform of real and even function results in a real function). To model the stellar photo- 
sphere we use a circularly symmetric uniform disk brightness distribution. The normalized 
visibility amplitude for a uniform disk can be written as: 



where Ji denotes the Bessel function of the first kind and first order, a is the angular diameter, 
and s is the radial spatial frequency (i.e., s = a/u^ + f^). Because the interferometric 
observations of P Cyg were obtained at sufficiently high spatial frequencies so that the Ha- 
emitting region was fully resolved, we can compare different models of brightness distribution 
for the envelope component. We choose to model our data with a uniform disk (UD) and 
a Gaussian distribution (GD), where in the former V^nv has the same functional form as 
equation E] except that a is replaced with 6'ud, and in the latter case 



where 6'gd corresponds to the full-width at half-maximum (FWHM) of the Gaussian model. 

The envelope component (Knv) in equation [H represented by either a uniform disk or a 
Gaussian radial distribution, was assumed to be circularly symmetric. This choice was based 
on a test performed on the observational data set being divided into 18 subsets based on 
the positional angles covered by the observations. Figure [6] shows the variation of the fitted 
diameter of a Gaussian model as a function of PA. Although a hint of variation is present in 
the data, more than 68% (based on itlo" expectation) of points fall within 5.5% of the mean 
value and thus we conclude that some of that variation might be solely due to the fact that 
the observations were acquired over many observing seasons and intrinsic variability at the 
5% level cannot be ruled out. Similar test performed on much smaller data set from only 




(2) 




(3) 
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one season does not reveal any variation beyond those expected based on the uncertainties 
associated with the fitted model parameters. Therefore, for the purpose of describing the 
geometrical structure of the Ha-emitting region of P Cyg, we decided to concentrate only 
on circularly symmetric models fitted to the V'^ data from the individual observing seasons, 
as well as to the entire data set. 

Table H] summarizes both UD and GD model results along with their corresponding 
reduced values (x^); which were used to assess the goodness-of-fit of each model. Based 
on the xt values obtained from model fits to the entire data set, it is clearly evident that 
the Gaussian model (with xl = 1-9) produced a significantly better fit to the data than a 
circularly symmetric uniform disk model (with xl = 2.8). However, neither the UD nor GD 
model fully reproduces the trend seen in the data obtained at high spatial frequencies as can 
be seen in Figured Although, the Gaussian model does produce a better fit, to fully account 
for the observational signature at the high spatial frequencies we require a model with an 
extra degree (s) of freedom. A natural extension to a Gaussian model is a two-component 
Gaussian model of the form 



2GD 



ki exp 



(vr^is) 



41n2 



+ (1 - ki) exp 



41n2 



(4) 



where 9i and 02 correspond to the FWHM values of the two Gaussian components, and ki is 
the fractional contribution fro m the first Gauss i an. S imila r two-componen t Gaus sian fits have 
been applied to other stars by lHofmann et al.l (120021 ) and lPreibisch et al.l (120031 ). Combining 
equations [H [2] and m and fitting the two-component Gaussian model to the observations from 
all three observing seasons results in angular diameters of 6i = 5.64 ± 0.17 mas and 62 = 
1.80 ± 0.11 mas, with a fractional contribution from the stellar photosphere, Cp = 0.70±0.01 
and ki = 0.57 ± 0.02 (i.e., with the two Gaussian components contributing approximately 
equally to the net Ha emission). The best-fit two-component Gaussian model had the 
lowest reduced value (yielding xl of 1.5 based on the entire data set) out of all the models 
considered in our analysis. The model curve corresponding to the two-component Gaussian 
model is also shown in Figure [21 

The success of the two-component Gaussian model to reproduce the observational sig- 
nature from the Ha channel is mostly due to the higher degree of freedom along the radial 
direction. In other words, one can in principle approximate any monotonically decreasing 
function as a sum of a large number of Gaussian functions. In fact, any other monotonically 
decreasing function can also be used instead of a Gaussian. For example, a sum of a Gaussian 
distribution and a uniform disk component of the form 



UD+GD 
env 



/i2 



Ji(7r6'uDs) 
7r6'uDS 



+ (1 - h) exp 



(7r6'GDg) 

41n2 



(5) 



- 8 - 



where ^ud and ^gd correspond to the angular diameters of the uniform disk and Gaussian 
components, respectively, and li is the fractional contribution from the uniform disk, yields 
equally good fit to the data as the two-component Gaussian model (resulting in a reduced 
of 1.5). The best-fit angular diameters are in this case: 6'ud = 3.06 ± 0.15 mas and 
6'gd = 5.46 ± 0.16 mas, with a fractional contribution from the stellar photosphere, Cp, of 
0.72 and h = 0.36. 



3.2. Central Star Diameter 



A central star diameter of 0.2 mas was assumed for our model (recall § 12. ip . Because 
varying the adopted diameter for the central star affects the calculated size of the Ha- 
emitting region, we test the dependency of the angular diameter of the Ha-emitting wind 
structure surrounding P Cyg on the angular diameter of the photospheric component. This 
is done by fitting equation [T] to the observations with the envelope component represented 
by a Gaussian model and the photospheric component with a UD model for a range of 
adopted stellar diameters. Using the observations from the 2007 observing season we find 
that the assumed central star diameter of P Cyg at most has only a 1% effect on the derived 
angular size of the wind structure when the adopted angular diameter of the central star is 
varied from to 1 mas. Therefore, an angular size of the continuum-emitting region that is 
underestimated by a factor of five results only in a 1 % effect on our best-fit angular diameter 
of the Ha-emitting region. 



3.3. Spectroscopic Results 



To test for any possible link between the changes in the strength of the Ha emission line 
and any variations in the size of the Ha-emitting region we examined the spectroscopic data 
for any variability on a time-scale comparable to our interferometric observations. We found 
that there was EW variability with an amplitude <10 A over a time scale of ^ 2 months, 
which is consistent with the short-term variability seen by lMarkova et al.l (1200161 ). However, 
the yearly mean Ha EW and yearly mean Ha peak signal showed no variation from year to 
year during our observing period suggesting relative stability in EW measure on a year-to- 
year time scale (see Table [3]). 

While our mean spectroscopic data suggest stability on a year-to-year scale, we cannot 
rule out variability of Ha em ission if such changes are correlated with changes in the con- 
tinuum level as suggested by iMarkova et al.l (l2001al ). If the continuum brightened and the 
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Ha has increased in strength by the same fractional amount, it would appear stable. Re- 
gardless of this apparent mean stability, we detect high frequency variations on time scales 
much shorter than one year. The upper panel in Figure O shows the Ha eq uivalent widths 



derived from our spectroscopic observations along with the data published by lMarkova et al. 



(12001 ah . 



Markova et al.l (12001 a ) obtained both spectroscopic and photometric observations of 
P Cyg and found that the EW measurements displayed a slow variability, with an amplitude 
of about 30 A and a duration of about 600 days, and a faster variability with an amplitude 
up to 10 A and duration of 40 to 60 days. While we were able to observe a variability of up to 
10 A in amplitude on a short (month-to-month) time scale, we did no t observe the Igng-terrn . 



30 A amplitude variability. If the long-term variability discussed in iMarkova et al.l ( 1200161 ) 
was present at the time of their observations, it is no longer present in our data. Therefore, 
based on our spectroscopic and interferometric observations, we conclude that P Cyg is in a 
quiescent phase of its lifetime, and it is not going through any drastic changes. 

The variability of the Ha EW in our data set is similar to the 10 A variation seen in 
the literature on a time scale of ~2 months. Based on our analysis of the equivalent width 
measures listed in Table [3l we conclude that the root mean square (RMS) variations are 
less than 10 %, and typical closer to the 2-5 % range. Therefore, for the purpose of our 
study we assume that the overall emission in the Ha line throughout our interferometric 
run was stable at the 10 % level. We also tested for variability in the Ha strength during a 
single night and found it to be virtually non- variable as well (see Table E]). The variations 
were within the expected observational uncertainties of ~3 %, mostly due to uncertainties 
associated with continuum level normalization. 



4. Discussion 

4.1. The Radial Intensity Falloff 

The deviation from a Gaussian shape of the spatially resolved wind structure of P Cyg 
is one of our most intriguing results. Examining the squared visibilities in Figure [U it 
is apparent that the observations at high spatial frequencies cannot be represented by a 
single- component Gaussian model. It is worthwhile to mention that the Gaussian intensity 
distribution is still a closer approximation to the radial intensity distribution from the wind 
structure than the uniform disk model (shown with dotted line in Fig. H]). To our knowledge, 
this is the first time the Ha-emitting region of P Cyg have been observed to clearly deviate 
from a Gaussian shape. 
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To reproduce our data at the highest spatial frequencies, we use the two-component 
Gaussian model as a way of parameterizing the spatially resolved signature of the Ha- 
emitting region of P Cyg. This suggests that the wind structure is indeed more complex in 
P Cyg as compared to Be stars. Because P Cyg po ssess radiatively d r iven s tellar wind, the 
velocity structure can be described by the /3-law of ICastor fc Lamerd (Il979l ). which implies 
that the wind is accelerated to its terminal velocity close to the star and then reaches a 
constant terminal velocity. This in turn means that the density structure of P Cyg will 
change at two different rates, one when the wind is still accelerating and one when the wind 
has already reached its terminal velocity at which point the density would be expected to 
drop with radial distance as r~^. Although this cannot be used as a direct expla nation for 
the t wo component structure because the Ha emission in P Cyg is optically thick (ILeitherer 
19881 ). by combining the effects of the density distribution with the radiative transfer effects it 
should be possible to use the data presented in this study to directly constrain a numerical 
wind model based on a /3 velocity law. Lastly, if the two-component wind structure is 
dominated by optically thick and thin effects, one could argue that the the inner (optically 
thick) region might be better represented by a uniform disk model and the outer region by a 
Gaussian model. This would be in agreement with the results presented in §3. II that showed 
that the observations can be equally well represented by a two-component Uniform-Disk and 
Gaussian model. 



P Cyg has been spatially resolved in the past using interferometric techniques. IVakili et al. 



(119971 ). using a single 17.7 m interferometric baseline, resolved the extended envelope of 
P Cyg at the Ha emission line. Although their data set was limited and did not allow them 
to fit a model more complex than a uniform disk model, the angular diameter of a uniform 
disk model they obtained was ^ud = 5.52 ± 0.47 mas. They also obtained complementary 
Ha spectrum on P Cyg on the same night they obtained interferometric observations with 
a Ha peak signal of 14.6 -FmaxZ-^cont- This is somewhat weaker than the ~20 -FmaxZ-^cont 
values we obtain based on our observations from 2005-2007 ( r ecall Fig. H]), although the 
effect of lower spectral resolution in the spectra of IVakili et al.l (119971 ) cannot be ruled out, 
which could expla in the lower peak-to-continuum ratio seen in the observation reported by 
Vakih et al.l Jl997[ ). 



Although a uniform disk model is co mpletely inconsiste nt with our observations (recall 
Fig. [2]), in order to compare our results to lVakili et al.l (119971 ) at effectively the same limiting 
spatial resolution, we fitted a uniform disk model to data from only shortest baselines (up to 
18.9 m). The best-fit uniform disk angular diameters ranged from 8.4 to 10.2 mas depending 
on which season was used in the fit (see Table Sj). Therefore, our UD diameter for the 
Ha-emitting wind structure is not consistent with the reported diame ter of 5.52 mas by 
Vakili et al.l (119971 ). This discrepancy might be due to the possibility that lVakili et al.l (119971 ) 
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has modeled the observational signal that contained both the mostly unresolved central star 
and Hof-emitting envelope with only one uniform disk diameter. This would be equivalent 
to our equation [H but with the Cp set to 0, which would tend to underestimate the angular 
extent of the Ha emitting region. The other possibility is that the Ha-emitting region 
actually grew in size fr om 1994 to 2 07. If we assume that the difference between our values 
and those reported by IVakili et al.l (119971 ) to be representative of a real change, then this 
results in an expansion rate of the Ha-emitting region (due to changes in opacity) over this 
thirteen year period of ~3 km s~^. This is well within the expected wind velocities of stars 
with radiation driven wind str uctures, especially P Cyg, which has a t erminal wind velocity 
almost two magnitudes larger (IBarlow et al.lll994j : iLamers et al.lll985l ). Furthermore, if the 
Ha-emitting region did indeed grow in size, we would expect the Ha emission to be stronger, 
which is supported by the fact that the -Fmax/-^c values in Table E] t hat range betweeii 17 
and 24 are larger than the ~ 15 peak value of the spectrum taken by IVakili et al. Jl997h . 



4.2. Signature of Asymmetry 

To test for the presence of a signature of deviation from point-symmetry in the signal 
in the Ha channel, we have plotted the weighted mean closure phases in Figure [31 If the Ha 
emission originates from an intensity distribution that is not symmetric across the origin of 
the photocenter, or equivalently the Ha emitting region is not concentric with the central 
star, we would expect to see a non-zero closure phase in the Ha channel. At a close inspection 
of Figure [3] we conclude that the Ha closure phases are generally very close to 0° . Even in 
the cases where the mean closure phase deviates from zero values, it is only at the level 
of a couple of degrees, which also happens to be only 2 to 3 times the uncertainty of the 
mean value. Therefore, we conclude that based on our observations we do not have a strong 
signature of deviation from point symmetry of the Ha-emitting region. 

Our conclusions are not necessarily inconsistent with the a very similar detection of 
a subtle ph a se va riation (at the level of ~ 30°) across the Ha emission line detected by 



Vakili et al.l (119971 ) who interpreted that as being produced by a localized spatial feature (i.e., 
a localized blob) within the structure of the wind. Assuming that the closure phase variations 
we detect across the Ha spectral channel in Figure [31 are caused by similar spatial feature 
(which is not necessarily expected since the observations were acquired at different epochs), 
we would indeed expect to see a much weaker signal across the Ha channel for two reasons. 
The first being that we measure the sums of three phases, which already have the tendency 
to lower the signal when negative and positive phases in a triangle are added together. The 
second effect is related to our much wider spectral channel that contains the entire Ha 
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emission line and a significant contribution from the central star, which most likely can be 
described well by a point-symmetric intensity distribution that contains only real components 
in the Fourier space (i.e., only real phases). The net result is that the real components will 
tend to lower the net complex phase detected in the Ha channel. 

Signatures of clumping in the circumstellar region of P Cyg on much larger spatial scales 
have also been reported. For exam ple, based on ne a r-diffr action limited observations using 
AO system on a 1.52 m telescope, IChesneau et al.l ( I2OOOI ) reported clumping at the scales 
of 200-600 mas. Similarly, P Cyg was also imaged using the Mu lti-Element Radio Linked 
Interferometer Network (MERLIN) array by lSkinner et al.l (119971 ). Their 6 cm observations 
produced images of the circumstellar structure on scales of 100-200 mas with a 50 mas 
resolution, reveali ng structural c hange s on a 40 day timescale along with flux variations 
at the 20% level. ISkinner et al.l (Il997l ) argued that the structural variations could not be 
attributed to variations in the mass-loss rate. They suggested recombination within the free- 
free emitting region responsible for the radio emission could explain their observations. The 
weak signature of the non-zero cl osure phase in th e NPOI dat a, if confirmed , could indicate 
that the dumpiness observed by IChesneau et al.l (I2OOOI ) and ISkinner et al.l ( 119971 ) already 
originates on scales of less than 10 mas (less than 25 stellar radii). 



4.3. Photometric Variability 

Lastly, to assess the level of photometric variability during the time frame covered by our 
interferometric observations, we used photoelectric V-band observations from the American 
Association of Variable Star Observers. The V-band light curve from 2005 to 2009 is shown 
in Figure where we only include data for which the measured check star (HD 193369) 
magnitude was no more than 0.04 mag away from the mean value of 5.573 mag. We see 
both a mean seasonal-timescale variability on the order of ~ 0.02 mag, and very short 
term variability (for observations taken within days to weeks) at the level of ±0.05 mag. 
Although the short-term variability appears to have large scatter, the photometric quality 
of the AAVSO PEP data is generally good, with internal errors on the order of 5-10 milli- 
magnitudes. Therefore, based on the photometric data we conclude that although the mean 
photometric level appears to change by less than 0.02 mag, short term variations upwards 
of 0.1 mag cannot be ruled out. This implies that the continuum level in P Cyg could vary 
up to 10 % level, which would affect directly the EW measures at the same level. This 
also strongly suggests that the variations in EW measures seen in Table |3] are caused by 
variations in the continuum level and not the emission component itself. 
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5. Summciry 

In general, mass-loss rates and stellar wind properties are extremely important charac- 
teristics of LBVs. Mass-loss rates for P Cyg have already been estimated through several 
methods, but none of them rely on spatially resolved Ha-emitting inner region. In this obser- 
vational study we have presented the most comprehensive interferometric data set on P Cyg 
to date. Our results provide a new and independent source of observational constraints that 
can be utilized when determining the properties of the wind structure in the context of the 
predicted Ha intensity distribution on the sky. 

We have fitted circularly symmetric uniform disk, Gaussian, and two-component mod- 
els to the interferometric observations and have shown how the Ha-emitting wind structure 
of P Cyg cannot be represented by a simple one component model. We conclude that the 
wind stnictnrc can be represented fully by a two-component Gaussian model with angular 
diameters of 9i = 5.64 ± 0.21 mas and 02 = 1.80 ± 0.13 mas. The interferometric sig- 
nature can also be equally well represented with a two-component model consisting of a 
uniform disk and a Gaussian intensity distribution. In that case the angular diameters are 
OvD = 3.06 ± 0.15 mas and Ogb = 5.46 ± 0.16 mas. While we cannot conclude if a double 
Gaussian model or a uniform-disk plus a Gaussian model better represents the wind struc- 
ture of P Cyg, it is clear that both models represent the wind structure better than a single 
Gaussian or uniform disk model. We conclude that P Cyg's wind structure is complex, pos- 
sibly containing multiple layers of varying opacities. The data presented in this study might 
be used to directly constrain wind models of P Cyg, perhaps allowing for an independent 
observational constraint on the /3-law. 

Based on our spectroscopic and interferometric observations we conclude that P Cyg 
is stable at the 10 % level. Based on the photometric observations we also rule out the 
possibility that large changes in the Ha EW measure were masked by changes in the contin- 
uum level. The lack of variability at a significant level combined with the spatially resolved 
Ha-emitting region that appears to be stable between the observing seasons is suggestive of 
an extra-quiescent phase of P Cyg between 2005 and 2008. 

The Navy Prototype Optical Interferometer is a joint project of the Naval Research Laboratory 
and the U.S. Naval Observatory, in cooperation with Lowell Observatory, and is funded by the 
Office of Naval Research and the Oceanographer of the Navy. We thank Anatoly Miroshnichenko 
and Glen Williams for valuable suggestions on how to improve the modeling section of this work. 
We are very grateful for the very detailed and constructive comments and suggestions we have 
received from the anonymous referee. We would also like to thank the Lowell Observatory for the 
generous telescope time allocation on the John S. Hall Telescope and Nevena Markova for providing 
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Corporation for Science Advancement. This research has made use of NASA's Astrophysics Data 
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Fig. 1. — Sampling of the (M,w)-plane in the Ha channel for P Cyg on multiple baselines 
ranging in length from 15.8 to 79.3 m obtained on thirty-seven nights of observation in 2005, 
2007 and 2008 {open circles). Possible coverage at each baseline is also shown from 6 hr east 
of meridian {dotted line) to 6 hr west of meridian {solid line). 
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Fig. 2. — The calibrated squared visibilities from the Ha channel of P Cyg obtained on 11 
unique baselines from 37 nights of observation in 2005, 2007 and 2008. A signature of a cen- 
tral star represented by a UD model with angular diameter of 0.2 mas is shown (dashed-line) 
along with three models representing the Ha-emitting envelope. The envelope component is 
modeled with a uniform disk {dotted-line) , a single component Gaussian {dash-dotted-line) , 
and a two-component Gaussian model {solid-line) . 
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Fig. 3. — Weighted mean closure phases (with a quaiiratic trend removed from the contin- 
uum channels) obtained at each possible baseline triangle obtained in 2005, 2007 and 2008 
observing seasons. The signal from the channel containing the B.a emission is marked with 
red squares. The closure phases containing the AC-AE-AN and E6-AE-AN stations in the 
2007 season are plotted in two panels since one of the baselines (AE-AN) was observed at 
two independent output beams from the beamcombiner (indicated in the square brackets 
with a SI or S2 designator), thereby allowing calculation of two closure phase quantities. 
The error bars are based on the uncertainty of the weighted average. 




Fig. 4. — A representative Ha line profile of P Cyg obtained on 2007 June 22. The equivalent 
width of the emission line is -80.0 A. 
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Fig. 5. — Uyyer panel: P C yg Ha equivalent widths from 1988 Jul to 2007 Dec with data 
from iMarkova et al.l ( 12001 d ) (triangles) and this study (squares) shown. The vertical lines 
mark the periods when interferometric observations were obtained at the N POI (dashed 



lines) , and the single night of interferometric observation of P Cyg reported by IVakili et al. 



( 119971 ) (dotted line). Lower panel: The AAVSO Johnson V-band photoelectric light curve of 
P Cyg for the period 2005-2009, obtained based on measurements relative to the comparison 
star HD 188892 (V = 4.936, B - V = -0.086). Only data where the measured check star 
magnitude (HD 193369; V = 5.573) lies within ±0.04 mag of the known value are plotted. 
Typical photometric errors per point are 5-10 milli-magnitudes. The smoothed moving 
average (with a 1000 d window) is also shown (solid line) to guide the eye for possible 
long-term variability. 
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Fig. 6. — Diameter of a Gaussian model fitted to subsets of squared visibility data divided 
based on ranges of PA on the sky. Each point represents a 10° range in PA. The mean 
diameter of 4.29±0.03 mas based on a model fit to all data {dashed line) and 5.5% variation 
around that value {dotted lines) are also shown. 
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Table 1: NPOI Observing 


; Log of P Cyg 


UT Date 


^ of Ho; squared 




visibilities 


2005 Jun 2 


48 


2005 Jun 3 


19 


2005 Jun 4 


48 


2005 Jun 6 


18 


2005 Jun 8 


20 


2005 Jun 9 


84 


2005 Jun 11 


19 


2005 Jun 13 


50 


2005 Jun 16 


90 


2007 Jul 1 


20 


2007 Jul 2 


20 


2007 Jul 4 


20 


2007 Jul 12 


15 


2007 Aue 8 


24 


2007 Aug 9 


16 


2007 Aug 10 


66 


2007 Aug 11 


56 


2007 Aug 12 


28 


2007 Aue 18 


72 


2007 Aug 19 


31 


2007 Aug 20 


82 


2007 Aug 21 


33 


2007 Aug 22 


92 


2007 Aug 23 


17 


2007 Aug 24 


12 


2007 Aug 25 


41 


2007 Aug 27 


12 


2007 Aug 29 


12 


2007 Sep 5 


36 


OfinS Tnn 97 


1 fi 


2008 Jun 28 


40 


2008 Jun 30 


36 


2008 Jul 2 


28 


2008 Jul 3 


60 


2008 Jul 7 


64 


2008 Jul 8 


89 


2008 Jul 9 


100 



Total 



1534 
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Table 2: Calibrated Ha Squared Visibilities of P Cyg 



Julian Date Spatial Frequency u Spatial Frequency v 



(JD - 2,450,000) 


(10^ cycles/radian) 


flO^ cycles /radian) 


Ho 


Baseline 




300 


-18.364 


19.654 


0.666 ± 0.028 


AC-AE 


3523.: 


863 


— 11. lO^ 




fi «9Q _|_ fi fiAD 
u.uzo HZ u.u^u 


AP AN 


3523.: 


863 


-7.180 


52.326 


0.615 ± 0.017 


AE-AN 


3523.: 


863 


-37.854 


29.298 


0.567 ± 0.023 


AC-E6 


3523.: 


863 


67.722 


-25.810 


0.512 ± 0.026 


AW-E6 


3523.: 


863 


29.868 


3.488 


0.671 ± 0.045 


AC-AW 


3523.: 


891 


-21.636 


17.482 


0.705 ± 0.023 


AC-AE 


3523.: 


891 


-8.106 


-33.720 


0.633 ± 0.039 


AC-AN 


3523.: 


891 


-13.530 


51.202 


0.515 ± 0.023 


AE-AN 


3523.: 


891 


-43.328 


24.890 


0.526 ± 0.025 


AC-E6 


3523.: 


891 


74.493 


-18.088 


0.501 ± 0.033 


AW-E6 



Note. — Table[2]is published in its entirety in the electronic edition of the Astronomical Journal. 



-25 - 



Tal)le 3: Spectre )sc()i)ic Yaria.l)ilit}' for P Cyg 



UT Date 


-EW (A) 


% Changet 


^max 
Fr. 


% Changet 


2005 Apr 1 


82.1 


3.5 


20.8 


6.7 


2005 May 19 


74.9 


5.7 


18.7 


4.0 


2005 Jun 28 


81.0 


2.0 


20.8 


6.7 


2005 Jul 16 


80.3 


1.1 


20.8 


6.9 


2005 Aug 22 


79.7 


0.5 


19.7 


1.1 


2005 Sep 16 


78.8 


0.8 


19.1 


2.0 


2005 Sep 17 


81.3 


2.4 


19.6 


0.6 


2005 Sep 27 


81.0 


2.1 


19.5 


0.1 


2005 Oct 11 


81.0 


2.1 


18.9 


3.1 


2005 Nov 8 


73.6 


7.3 


17.0 


12.9 


2005 Mean 


79.4 ± 0.9 




19.5 ± 0.4 




2006 Jun 10 


81.4 


7.3 


19.7 


7.3 


2006 Aug 3 


70.1 


7.7 


16.6 


9.4 


2006 Sep 6 


75.3 


0.8 


17.9 


2.7 


2006 Oct 3 


70.3 


7.3 


17.2 


6.2 


2006 Oct 28 


78.1 


2.9 


19.0 


3.3 


2006 Oct 31 


77.4 


2.0 


19.0 


3.4 


2006 Nov 1 


78.6 


3.5 


19.2 


4.4 


2006 Mean 


75.9 ± 1.6 




18.4 ± 0.4 




2007 May 3 


70.5 


10.8 


17.3 


17.7 


2007 Jun 4 


75.8 


4.1 


18.7 


11.1 


2007 Jun 22 


80.0 


1.2 


19.7 


6.4 


2007 Oct 23 


83.4 


5.5 


22.5 


6.8 


2007 Oct 24 


78.0 


1.3 


21.1 


0.3 


2007 Oct 27 


76.5 


3.2 


20.9 


0.6 


2007 Oct 28 


74.7 


5.6 


20.2 


3.9 


2007 Dec 18 


80.9 


2.4 


22.6 


7.5 


2007 Dec 19 


85.5 


8.2 


23.7 


12.7 


2007 Dec 20 


85.2 


7.8 


23.7 


12.5 


2007 Mean 


79.0 ± 1.6 




21.0 ± 0.7 





^ Changes calculated with respect to yearly mean EW and peak strength. 
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Table 4: Best-Fit Model Parameters for P Cyg 



Model Season Diameter (mas) 2°*^ Diameter (mas) fci or /i Cp xt 

Uniform Diskt 2005 9.41 ± 2.50 77. 77. 0.84 ± 0.03 

Uniform Disk 2005 6.74 ± 0.05 ... ... 0.78 ± 0.01 2.6 

Gaussian 2005 3.97 ± 0.05 ... ... 0.75 ± 0.01 1.5 

Double-Gaussian 2005 8.7 ± 5.3 2.08 ± 0.15 0.41 ± 0.03 0.70 ± 0.01 1.2 

Uniform Diskt 2007 10.17 ± 0.32 77. 0.84 ± 0.01 

Uniform Disk 2007 7.37 ± 0.03 ... ... 0.79 ± 0.01 3.7 

Gaussian 2007 4.37 ± 0.03 ... ... 0.76 ± 0.01 2.4 

Double-Gaussian 2007 5.33 ± 0.15 1.53 ± 0.15 0.60 ± 0.02 0.70 ± 0.01 1.7 

Uniform Diskt 2008 8.40 ± 1.50 77. 77. 0.79 ± 0.04 

Uniform Disk 2008 7.06 ± 0.06 ... ... 0.78 ± 0.01 1.5 

Gaussian 2008 4.37 ± 0.07 ... ... 0.75 ± 0.01 1.3 

Double-Gaussian* 2008 10 ± 9 2.17 ± 0.07 0.50 ± 0.02 0.70 1.1 

Uniform Diskt Ml 9.83 ± 0.26 77. 77. 0.84 ± 0.01 

Uniform Disk All 7.21 ± 0.03 ... ... 0.79 ± 0.01 2.8 

Gaussian All 4.29 ± 0.03 ... ... 0.76 ± 0.01 1.9 

Double-Gaussian All 5.64 ± 0.17 1.80 ± 0.11 0.57 ± 0.02 0.70 ± 0.01 1.5 



Uniform Disk + Gaussian All 3.06 ± 0.15 5.46 ± 0.16 0.36 ± 0.02 0.72 ± 0.01 1.5 

t Models fitted to data for baselines up to 18.9 m. 

i There were insufScient data at the longest baselines to fit for four model parameters and instead a fit was obtained 
with Cp fixed at 0.70. 
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Table 5: Single Niglit Spcctrc^copic Variability' for P Cyg 



UT Date 


-EW (A) 


% Change^ 


max 


% Change^ 


2007 Dec 20 


85.2 


0.1 


23.7 


0.0 


2007 Dec 20 


85.9 


0.9 


23.8 


0.7 


2007 Dec 20 


84.4 


0.9 


23.5 


0.7 


2007 Dec 20 


85.6 


0.5 


23.8 


0.4 


2007 Dec 20 


85.2 


0.1 


23.7 


0.0 


2007 Dec 20 


85.8 


0.7 


23.8 


0.6 


2007 Dec 20 


86.0 


1.0 


23.8 


0.6 


2007 Dec 20 


83.4 


2.1 


23.4 


1.2 


2007 Dec 20 


85.0 


0.2 


23.6 


0.2 


2007 Dec 20 


85.2 


0.0 


23.7 


0.0 


2007 Dec 20 


83.4 


2.1 


23.4 


1.2 


2007 Dec 20 


84.7 


0.6 


23.6 


0.4 


2007 Dec 20 


85.3 


0.1 


23.6 


0.1 


2007 Dec 20 


85.5 


0.4 


23.7 


0.3 


2007 Dec 20 


85.4 


0.2 


23.7 


0.1 


2007 Dec 20 


85.7 


0.6 


23.8 


0.5 


2007 Dec 20 


86.1 


1.1 


23.8 


0.6 


2007 Dec 20 


85.5 


0.4 


23.7 


0.3 


2007 Dec 20 


85.1 


0.1 


23.6 


0.3 


2007 Dec 20 Mean 


85.2 ± 0.2 




23.67 ± 0.03 




Changes calculated with respect to 


nightly mean. 


All spectra were 


obtained with 



180 s of integration. 



